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We consider the orbitally degenerate 3-band Hubbard model with on-site interactions which favor 
low spin and low orbital angular momentum using standard second order perturbation theory in the 
large Hubbard-U limit. At even integer filling this model is a Mott insulator with a non-degenerate 
ground state that allows for a simple description of particle-hole excitations as well as gapped spin 
and orbital modes. We find that the Mott gap is generally indirect and that the single particle 
spectrum at low doping reappears close to even filling but rescaled by a factor 2/3 or 1/3. The 
model captures the basic phenomenology of the Mott insulating and metallic fuUerides A^Ceo- This 
includes the existence of a smaller spin gap and larger charge gap at even integer filling, the fact 
that odd integer stoichiometries are generally metallic while even are insulating, as well as the rapid 
suppression of the density of states and superconducting transition temperatures with doping away 
from a; = 3. 

PACS numbers: 



I. INTRODUCTION 

The physics of Mott insulators has emerged as a key ingredient in the study of strongly correlated systemai, largely 
motivated by the understanding that much of the exotic physics of the underdoped cuprate superconductors has 
connections to the undoped Mott insulating state. The basic model for a Mott insulator discussed in connection 
with the cuprates is the large-U half-filled Hubbard model. For this model the low energy spectrum at half-filling is 
reasonably well understood, but the higher energy spectrum as well as the physics away from half-filling is still very 
much an open problem. Much of the difficulty seems to trace back to the fact that the Mott insulating ground state 
for this system is not known. 

Here we will study a model of a Mott insulator which has a simple non-degenerate ground state and where we can 
describe the spectrum of excited states and explore the physics at and near the Mott transition in a well controlled 
manner. The model is a degenerate three channel Hubbard model which has multiplet splitting on-site interactions 
that favor low spin and low orbital angular momentum. This model is the simplest model with spin and orbital 
symmetry which at even integer filling allows for a Mott insulating non-degenerate ground stated Since the ground 
state is non-magnetic, the spin physics which complicates the single band Hubbard model is absent at low energies. 
Although interesting in itself as a natural and simplifying extension of the Hubbard model, our main motivation for 
studying this model is that it very naturally captures some of the most distinctive phenomenology of the alkali doped 
Ceo compounds, the fuUerides. 

Crystalline Ceo is a band insulator with a completely filled molecular orbital^ Doping with alkali atoms, forming 
AnCgQ, transfers electrons into the lowest unoccupied molecular orbital (LUMO) which according to elementary 
molecular theory is threefold degenerate. These form three spin degenerate bands according to the band theory of 
the crystal. It is believed that the charge transfer is complete and that the influence of the alkali ions is in general 
negligible apart from changing the crystal structure and the corresponding band theory. 

If these were simple metals, band theory would would predict a metallic state for any doping < rt < 6 and 
possibly superconductivity with transition temperatures which would follow roughly the density of states at the Fermi 
energy as the doping is varied. Experimentally, however, it turns out that the compounds with even integer doping 
(n = 2,4) are non-magnetic insulators^i^ and superconductivity is only seen in a narrow range around half-filling 
(n = 3) and with transition temperatures that are sharply maximized here?. Band theory augmented by BCS theory 
of superconductivity fails to reproduce this behavior. 

It is well known that correlation effects are important as these are weakly bound molecular solids with a narrow band 
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width of approximately .5eV and it is believed that the on-site Coulomb repulsion is approximately leV. One might 
therefore expect that the system would be a Mott insulator at any integer filling. This, however, does not explain 
why systems with even integer filling are generally insulating while systems with odd integer filling are generally 
metallic. A second related question is why the insulating materials at even integer filling are non-magnetic. From 
Hund's rule, we would expect the highest spin configuration to be the molecular ground state, S—1 for n=2 and n— 4, 
and consequently some sort of magnetic ground state is to be expected for the solid. 

One explanation for the violation of Hund's rule in the insulating materials is that the Jahn- Teller (JT) effect 
counteracts Hund's rule, with the lowest energy JT distorted state of the Cgo molecule at even integer occupation 
being a spin singlet.- However, a problem with the naive Jahn- Teller scenario is that static distortions of the Ceo 
molecules have not been detected in solid Cgcn. Indirect evidence has come from the existence of two energy gap 
scales in the insulating systems, a smaller spin gap of around bOmeV and a larger charge gap of the around bOOmeV. 
The large charge gap is quite clearly a Mott gap related to the intramolecular electron-electron repulsion, whereas 
the smaller spin gap has been linked to a singlet to triplet gap of JT distorted molecules,'^. 

A different scenario for the violation of Hund's rule suggested by Chakravarty et aii^^ and Baskaran and Tosattiii 
claims that electronic correlations on a Cgo molecule can break the degeneracy of any partially filled molecular orbital 
in such a way as to minimize the spin and orbital angular momentum. Chakravarty and coworkers looked at a Hubbard 
model on single Cgo molecule, with strong electron-electron repulsion (U) on each carbon atom, using second order 
perturbation theory in U. The linear in U terms give Hund's rule with a preferred highest spin configuration which 
minimizes the overlap between electrons. The second order term on the other hand prefers low spin configurations 
with large overlaps which can take better advantage of virtual excitations of the core electronSfi^ The validity of the 
qualitative features of the perturbation theory is supported by exact diagonalization of smaller Hubbard clusters^ 

One important feature of the scenario based on intramolecular electronic correlations is that in contrast to the static 
Jahn- Teller distortions the orbital symmetry of the molecule is preserved. More recently, it has also been emphasized 
in a series of papers by Tosatti and coworkersi^ti^JiSiii that Jahn- Teller phonons treated in the anti-adiabatic limit 
will also give rise to an effective electronic Hamiltonian with inverse Hund's rule which preserves the orbital symmetry 
of the molecule. It may thus be difficult to distinguish the latter scenario from the one based on electronic correlations 
and a full treatment of both the electron-electron and the electron-lattice interactions to decide which energy scales 
dominate is asked for. Along these lines a recent density-functional calculation docs find that Hund's rule is valid 
for an isolated Cgo and only counteracted by the JT effect^^, although issues such as electronic screening by the 
surrounding molecules may also be important siS. 

Regardless of the mechanism behind the inversion of Hund's rule on the Ceo molecule it is an interesting problem 
to study the effective model of solid Ceo which this naturally gives rise to, namely the three channel Hubbard model 
with "anti-Hund's rule" interactions. Here we study this model in the limit of strong interactions where we can do 
standard second order perturbation theory in the intermolecular hopping. We find that the even/odd integer doping 
effect, the existence of two energy gaps in the insulating systems, as well as the rapid variation of Tc with doping 
away from n — 3 are all natural consequences of this model. 

We will focus primarily on even integer filling (n=2 or 4) where the problem simplifies significantly because of 
the resultant non-degenerate strong coupling ground state. Based on this we can describe the spectrum consisting 
in general of an indirect Mott gap to particle-hole excitations and chargeless spin and orbital modes with a distinct 
energy gaps. We also derive the spectrum of a single particle or hole doped into the non-magnetic Mott insulator. 
The particle and hole spectra turns out to be exactly the same as the non-interacting band structure up to an overall 
rescaling of the hopping by a factor 1/3 or 2/3. These characteristic properties of the non- magnetic Mott insulator 
are in agreement with previous results on the same mode h^^i^^ 

Apart from the correspondence with the physics of the Mott insulating A2Ceo and A4C60 our main conclusion is 
that Ceo with a filling n = 2-|-xor4-|-a;of the LUMO band, with |a;| << 1, is best understood as a doped Mott 
insulator, manifested by a charge carrier concentration given by |x| instead of n. In addition the density of states is 
up to a rescaling factor the same as that of the non-interacting band structure at the band edges n=|a;|orn = 6— |a;| 
for particle and hole doping respectively. If we heuristically extend these results to larger a; « 1 it follows naturally 
the observed variation of Tc with doping for 2 < n < A within any weak coupling BCS like scenario as a consequence 
of the rapid drop in the density of states close to non-interacting band edges. 

We will not, however, discuss the actual mechanism of superconductivity, i.e. the source of pairing. In particular, 
we will not assume that the intramolecular singlet triplet gap is necessarily large enough to overcome the Coulomb 
repulsion and give rise to a bare attraction along the lines of the earlier workiiSiii In fact, a fit of the parameters of 
the model to experiment suggests that this is not the case. Nevertheless, any other mechanism, such as one based on 
attraction from electron-phonons interactions, certainly needs to include the strong electron correlations which are 
present as indicated most clearly by existence of Mott insulating phases. 

The paper is organized as follows: In Sec. ^Iwe define the model and the strong coupling limit we will primarily 
consider. Then we derive a perturbative effective Hamiltonian valid in the strong coupling limit. In Sec. Illll we study 
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the Mott insulator at even integer filling and describe the spectrum of excited states. We also discuss the distinction 
between even and odd integer doping and why at odd integer doping the system is more likely to be on the metallic 
side of the Mott transition. Then we use the results derived for the model to get an estimate of the parameters by 
comparing to experiment on A2C60 and A4C60. In Sec. IIVI we study the doped Mott insulator at a filling close to 
even integer and speculate on the implications of these results to the whole doping range 2 < n < 4 as well as the 
properties at higher temperatures. Finally, in Sec. 0we conclude. 



II. THE MODEL 



We will be studying a three band Hubbard model. The electrons occupy three degenerate p-orbitals {L — 1) on 
every site of the lattice and we define electron creation and destruction operators cl 
quantum number and spin respectively. We will be working in the basis where / 

The Haniiltonian reads 



^ and Cj, with site index, orbital 
-1,0', 1. 



with a nearest neighbor hopping 



and on-site interaction 



H = h + Hi 



{rr') 



(1) 



(2) 



Hi = 2^ -Unf. + JlK + JsS^ 



(3) 



with U, Jl-.Js > 0. Here 



Lr — J2ws^rAs^l,l-'^r,l's 



(4) 
(5) 
(6) 



are the number, orbital angular momentum and spin operators respectively. 

The lattice is three dimensional and we assume a point group symmetry which is such that the hopping preserves 
the threefold orbital degeneracy, although this may be relaxed as long as the symmetry breaking crystal field is small 
compared to the interactions. We are going to be studying this in the strong coupling limit defined as 



U » Js » Jl »t, 



(7) 



where t = max(t[;7 ). For Ceo band structure calculationsSS give t ~ lOOmeV and values of U of the order leV have 
been suggested^. A subsequent fit of the parameters of the model to experiment on AnCgo suggests that the large U 
limit is quite well satisfied with U/t ^ 5, while Js ^ Jl ^ t. The limit U » t is essential for our treatment while 
relaxing the other limits probably will not change the qualitative features of our results. 

At t = we find that the Hamiltonian is simply diagonalized in terms of the states \n, L, S, , S^)r with n 
the number of electrons on a site and (L, S) the total angular momentum and spin and the energy E{n^ L, S) = 
\Un'^ + JlL{L + 1) + JsS{S + 1). In the strong couphng limit the effect of the hopping t will be to reduce the 
translational degeneracy It is the topic of the next section to derive an effective Hamiltonian which describes these 
low energy degrees of freedom. 



A. Effective strong coupling Hamiltonian. 

Let us define the short-hand notation p = {n,L,S) and a = {L^,S^). As shown in Tabled any single site states 
\p,a)r and \p',a')r have the same energy, barring accidental degeneracies, if and only Hp — p' . This implies that 
any eigenstate of Hj will not mix different representations at a site and can be written as Or V'p qIPi Q;)r, where 
V'p a ^ only for one particular p{r). Consider now the action of the hopping term h on an eigenstate. This will 

affect two arbitrary nearest neighbor sites r and r', giving \p, a) r\p' ,€('),.' ^ '^qBq'B' P)rW , P')r' , where M are 
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(0,0) 
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(1.5) 


5(7 + 2Jl + |Js 
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(0,0), (2,0), (1, 1) 


2(7 + {0, 6Jl, 2Jl + 2Js\ 
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(l,i),(2,i),(0,f) 


lU + {2Jl + f Js, 6Jt + f Js, Js} 


4 


(0,0), (2,0), (1,1) 


8(7 + {0, 6Jl, 2Ji + 2Js} 


5 




f (7+ fJs 


6 


(0,0) 


18U 



TABLE I: Spectrum of the interaction, Eq. |21 at a single site with states specified by occupation n, total orbital angular 
momentum L and total spin S. 



the matrix elements of h between the in and out states. Unless the initial and final nearest neighbor states have 
only exchanged particle number and representation between the sites, i.e. q = p' and q' — p, there will be a non- 
zero energy difference AE = E{q) + E{q') — E{p) — E{p') with respect to Hj which is some linear combination of 
U,Jl,Js- In the strong coupling limit Eq. there is no accidental degeneracy and the energy difference satisfies 
AE >> t, which will allow us to do a perturbation theory in t/ AE. Henceforth we will use the AE to define the 
absolute value of the minimum non-zero energy difference between two different nearest neighbor state configurations. 
From table ^ we can check that in the strong coupling limit this corresponds to AE = 4:Jl, given for instance by 
E{3, 2, i) + E{2, 0, 0) - Ei2, 0, 0) - E{3, 1, i). 

We will be using the Effective Hamiltonian method described as applied to the one band Hubbard model in, for 
instance, FazekasSi. We introduce a canonical transformation 

j2 

Heff = e'^He^'^ + i[S,H] + -[5, [5, H]] + ... , (8) 

where the goal is to find "Heff such that it does not couple the different energy sectors of Hj. 

Let us consider h as a matrix which acts on nearest neighbor states \p, a)r\p' , a')r' = \a)rr' which we in a short-hand 
notation will denote by a single index a,b,c, .... We write 

/i= ^/.-'^ ^/,-'|a)„,(&|„,, (9) 

{rr') {rr') 

where 

Kl = Hrr'Y.twicl^lsCr',l's + h.C.)|6)„. . (10) 

w 

Next we split h^l into h^l = h'^'^^ + h]^^^ where h^^^^ — ^E(a).E{b), only connects states with the same energy 

and hl^^^ = h^l — h^^^ which connects states with different energy. It follows that we can also write h = h'^ + h^. 
We introduce S = Si + S2 where Si - 0{t/AE) and ^2 ~ 0{t^/AE^) and expand Eq. lito 0{t^/AE), 

Heff = Hi + h'> + h'+ i[Si,Hi] +i[Si,h° + h'] + '^[Si, [Si, Hi]] + i[S2,Hi] . (11) 

To first order in t we want to cancel by taking Si to solve 

i[Si,Hi] = -h^ (12) 

Clearly Si needs to act only on nearest neighbor sites. Defining Si — J2{rr'}^i '' ^ 12 {rr') ^i^abl'^) rr' {b[rr' gives 
from Equation H12|l 

^J2,,l,rr' ^ J2 [Sr',Hf] 

{rr') (^j^^f'^j^" 

= Y.SZj\a)rr'{b[rr',H^I+H^'] 

{rr') 

= J2^Zabim^E{a))\a)rr'{b[rr', (13) 

{rr') 
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where we have used the fact that we work in the eigenbasis of Hj such that (iJJ + Hj )\a)rr' = E{a)\a)rr'- We thus 
arrive at the solution 

^ E{b)^-E{a) ■ ^^^^ 

Recall that h]^^^ is defined as to only have non-zero entries when \E{b) — E{a)\ > AE so that Si ^ t/ AE. 
Given the solution for Si obeying Eq. E|we can rewrite Eq. 1111 as 

Heff =Hi+ h" + i[Si,h^^ + h^] - f h^] + i[S2, Hj] 

= Hi + h'' + /jO] + ^[Si,h^] + i[S2,Hi] (15) 

In analogy to using Si to cancel terms linear in t which connect different energy subsectors we can now use 52 to 
cancel similar terms to order which arise from the commutators of Si and h in Eq. E| Technically this is slightly 
more involved because there is also three site next nearest neighbor interactions generated to second order in t and 
we refer the reader to the Appendix for the details. 

The final expression for the effective strong coupling Hamiltonian reads 

Heff = Hi + J2 i^effY:b\a)rr'{b\rr'+ ^ (Hef fVat"-" \a),r'r" {b\rr'r" + 0(1^ AE^) , (16) 

[r.r') ,ab {r,r'r") ,ab 

with 

Kl - E(cT- E(a) I '^-^''^ 

(nj \rr'r" ^ '^ac '^cb '^ac '^cb r /1 o^ 

(K//)., 2. Eic)-Eia) ^^(-^''^W 

c,E{c)=^E{a) ^ ' ^ ' 

Here \a)rr' = {n, L, S, , S^)r\n' , L' , S' , V , S'^)r' and \a)rr'r" are arbitrary nearest neighbor and next nearest neigh- 
bor eigenstates of Hj with energy E{a) and h^l is defined according to Eq. EH with a straightforward extension for 
three site states. 



III. THE MOTT INSULATOR 



In the strong coupling limit, Eq. |7| of this model the system is an insulator at any integer filling with a gap to charge 
carrying excitations of order U. At odd integer filling the ground state is the highly degenerate \GS, 1/3/5) = Ylr = 
1/3/5, L = 1, 5 = i, L^, 5"^)^ and the hopping t will introduce strong correlations in analogy with the half-filled single 
band Hubbard model. We will return to this problem briefly in discussing the fact that A3C60 is generally on the 
metallic side of the Mott transition. 

However, at even integer filling, n = 2 or 4, the ground state is non-degenerate, \GS,2/A) — Hr I'* ~ 2/4, L — 
0, S = 0)r- Here, we will find that much of the physics can be understood in terms of a simple non-interacting single 
particle picture. This will be the focus of the subsequent discussion and the main issue of the paper. We will discuss 
n = 4 to compare with experiments on K4C60 and Rb4C6o but the treatment of n = 2 is completely analogous. 

At filling n = 4 there are three types of excitations of the ground state that can be readily identified and 
which are indicated in figure ^ We can excite a single site into a higher energy multiplet, creating states 
|4,2,0,L^),nr'^r|4,0,0),. or |4, 1, 1, L^ 5^), Il.'^r K, 0, 0),. with energy A = E{4,2,0) - E{i,0,0) = 6Jl and 

A = -E(4, 1, 1) — E{4:, 0, 0) — 2Jl + 2Js respectively. There is also "particle-hole" excitations in different multiplets, 
the one with lowest energy being |3, 1, i, L^, S'^)r|5, 1, ^,L'\S"')r' Ilr"^r.r' 0, 0)^" with energy Aq = U+Ul + ^Js- 
The energy Aq is the Mott gap to O'th order in t. 

The degeneracy of the excited states will be lifted by the hopping h and using the effective Hamiltonian Eq. 1171 we 
can study the spectrum perturbatively in t. 



A. Spin gap, spin and orbital modes 

The spin and orbital excitations |4, 1, 1, L^, S^)r HrV 1^' ^' ^^^^ 1^' ^' ^^'^^ Ylr'=^r Oi 0)r' are degenerate in 
and as well as position r. Acting with the effective Hamiltonian Eq. 1171 will split the degeneracy in space (as 
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(n=3,L=l,S=l/2) + (n=5,L=l,S=l/2) 
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(n=4,L=l,S=l) 
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(n=4,L=2,S=0) 
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, ^0 (n=4,L=0,S=0) 



FIG. 1: Spectrum ai t = and doping n = 4 showing the lowest energy excitations of the three distinct kinds discussed in the 
text. 



well as the orbital degenercy) leading to a band description of these states. Clearly to first order in h (t), the effective 
Hamiltonian will not affect these states as it necessarily creates a high energy particle-hole state. For the same reason 
the second order three site interaction will not contribute. However, the second order nearest neighbor term can hop 
the excited state at r to a nearest neighbor r' or to the same site r through an intermediate particle-hole state. 
Focusing on the spinless excitations we define the "bosonic" operator 

= |4,2,0,m),(4,0,0|, (19) 

in terms of which we can write an excited orbital state as bl j^\GS, 4). These operators do not obey proper commutation 
relations but if we consider only single particle physics this is irrelevant. With these operators we can write the 
following single particle Hamiltonian describing the dynamics of such excitations 

H — ^ ^ ^mm'^r,m^r' ,7n' + ^ ^^ 'tmm'b\„^br^m' i (20) 
{r,r') r 

where 

^ ~U + 3/ij^'^^ (4,2,0,m|,(4,0,0|..t[;'(c;,,c,, ^,,+h.c.)|ph) x 

(ph|i-;(ct_^.^c,,^^.,, +h.c.)|4,0,0),|4,2,0,m'),, , (21) 

and similarly for tmm' ■ Here |ph) are particle-hole states |5, 1, ^,l,s)r\3, 1, ^,l',s)r' (or with reversed positions) and 
the sum is over all indices except m, m' and r, r' . Similar contributions but with larger denominators from intermediate 
particle hole excitations in higher energy multiplets have been ignored in Eq. [Til 

Deriving the Hamiltonian Eq. I2UI is thus a straightforward problem. The Hamiltonian is non-interacting and can in 
principle be diagonalizcd in momentum space. In general the five orbital components will be mixed with a complicated 
band structure. The same calculation can be carried out for the spinful (n — A, L — 1, S — 1) states, where there are 
nine coupled states per site. Although the band width of these spin and orbital modes are naively order the 
high degeneracy of these states is expected to give a significant broadening of band width, perhaps by a factor of the 
spin and orbital degeneracy of these states. The spinful excitations should be readily detectable experimentally as a 
band of magnons with some gap A^, a spin gap. We will return to the issue of the spin modes in discussing NMR on 
AjCeo and A4C6o- 
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B. Charge Excitations 

At filling n = 4 the lowest energy charge carrying excitations are the "particle-hole" states 
|3, 1, Z, s)r|5, 1, i J', s')r' nr"5^r r' 0, 0)^" with an energy Aq = U + AJ^ + fJg. Acting with the effective Hamil- 
tonian Eq. El on such a state will translate either the particle (n = 5) or the hole (n = 3) to a nearest neighbor site 
through the linear term in h (t). To second order, there is next nearest neighbor hopping through an intermediate 
state in a higher energy multiplet. Since it is higher order, 0{t'^ /AE), we will ignore this contribution. However, it 
should be noted that the denominator AE for this process is proportional to Js and/or Jl not U, implying of course 
that we should consider the higher energy particle and hole multiplets if the strong coupling limit Eq. [7|is not strictly 
valid. Note that due to the energy constraint in the effective Hamiltonian the particle and hole are not allowed to 
annihilate by hopping to the same site. We will ignore the constraint which in three dimensions is expected to give a 
vanishing contribution to the spectrum at low densities. 

1. Charged single particle excitations 

In order to understand charge transport, we need to consider the single particle and single hole states 
|5,Z,s)rnr' \^)r'^r and |3,Z,s)rnr' \'^)r'^r- Here we have dropped the multiplet indices (^,5") — (1, 5) or (0,0). 
When acted on by Tie/ / these excitations are translated to a nearest neighbor through h with a matrix element 

{4\r{iJ\s%,{Y,t7'A.J,<rCr',,',.)\3J,s)r\4)r' = -^^..'i^', (22) 

and similarly with a factor —1/3 for holes and 2/3 for particles with respect to the n — 2 ground state. The matrix 
elements in Eq. 1221 follow from the explicit expressions for the states 

15,1,^,^,5) = \/3cL|4,0,0) 

13,1,^,^,5) = y|cL|2,0,0)^y|(2.)^l'|-ic_,_j4,0,0) 

|l,l,i Z,s) = c|jO) = V3(25)2|'l-ic_;_j2,0,0), (23) 

which is a simple exercise in elementary quantum mechanics to derive. 

We note the reduced magnitude of the matrix elements compared to the single particle particle hopping on empty 
sites 

(0|,.(l,r,s|,„(^ ty;cl^^,^^Cr,,,^)\lJ,s')r\0)r' = <S..'tP^' . (24) 

This is a general result of the reduced Hilbert space due to constraining the n particle states to the lowest energy 
multiplet. 

Let us define the particle and hole creation operators 

4.rls = |5,Z,s)r(4|. 

clris = |3,Z,s).(4|., (25) 

through which the particle and hole states can be written as Cg ^.^JCS, 4) and Cg ^.^JCS, 4) with \GS, 4) = Yir |4, 0, >r- 
Note that these operators do not obey on-site anticommutation relations, we can only use them with confidence in 
describing non-interacting single-particle physics. 

In terms of these particle and hole operators we now straightforwardly arrive at the following single particle Hamil- 
tonians for the particle and hole states 

^5= E ^*n''4,WsC5,r'Z'.+E(9/2t^ + 3/4^S + 2JL-M)4,H.C5,H. (26) 
{rr'),ll's r,ls 

H3^- E fu'clrls^3.r'l's + Y.^-l/'^U + 3/4Js + 2 + /i)4,,;,C3,H. , (27) 
{rr'),ll's r,ls 
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where the on-site energy is defined with respect to the E{4, 0, 0) = 8f/ — and where we have introduced the chemical 
potential /i. 

These are just simple tight binding Hamiltonians which we can diagonalize in momentum space in terms of states 

\5,kJs)=cl,JGS,A) ^J^J2^'''%ris\GSA) 

V ^ 

\3,k,ls) = 4,,,JG5,4) ^j^Y.e-^^'-%^JGS). (28) 

In general these are not eigenstates due to the fact that the hopping is not diagonal in I and we would get three (or 
more if there are inequivalent sites) non-degenerate bands depending on the precise nature of the hopping integrals 
which depend on the crystal symmetry. Quite intriguingly, as seen from Equations 1261 and 1271 the band structure 
for these particle and hole excitations from the Mott insulating ground state is precisely the noninteracting band 
structure up to a rescaling factor. The reason for this remarkably simple behavior is that the ground state at even 
integer filling is in the trivial (L = 0, 5* = 0) representations of spin and angular momentum and as such is basically 
equivalent to the zero-particle vacuum. 



2. Mott gap 

The Mott gap is defined as the gap between the ground state and the lowest energy charge carrying excitation of 
the insulator. What is often measured however is the optical gap as defined by optical conductivity or reflectivity 
measurements and we will be interested in calculating this too. 

The Kubo formula for the optical conductivity which is the short wavelength limit of the electrical conductivity is 
at zero temperature 

where Im) are excited states with HuJm — E„i — Eqs-^'^ The current operator can be derived from the continuity 
equation i'h~^[n{r), H] -f V • j{r). For the Hamiltonian Eq. ^ considered here only the tight binding part contributes 
and we get 

poll's \ ^ 5 / 

where 5 is the set of nearest neighbor lattice vectors. For a cubic or orthorhombic lattice this simplifies to 

^^'i) = E ^^'<p- ^) i+,,is^P,i's ■ (31) 

p,ll' s,S 

The current operator creates particle- hole pairs when acting on the ground state 105,4) and we expect to get 
non-zero matrix elements with states 

|fc/s,fc'/V)p.h-4,,,4,,,,,,|G5,4) (32) 

defined according to Eq. 1281 and 

With this we can calculate the matrix element 

(fcZs,fcrs%.hr(g = 0)|G5) = -||(2sO'l''l"'^5"sin(^.jX-z''5fe,fe'^s,-.' (33) 

6 

Introducing an explicit band structure, i.e. defining t^f , we can in principle calculate the optical conductivity due 
to the particle hole excitations. A lower bound to the support of the sum in Eq. [5^1 will tell us the optical gap, below 
which cj, (7{lu) will decay. We simply maximize the kinetic energy of the particle and hole under the zero momentum 
constraint in the usual manner to obtain this. We will ignore the possible complications due to the dispersion and 
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angular momentum part 5° sin(fc • which may kill the matrix elements at some high symmetry points. 

Moving slightly away from such symmetry points will give a finite contribution to the response. 
The optical gap is thus given by 

Aoptical >U + AJl + + - ^4(^))min(£,.,,) ' (34) 

where e* (fc) is the kinetic energy of the j'th band at momentum k and min(/c, i,j) means minimizing with respect to 
the momentum and the band indices. 

To get the Mott gap we just need to find the lowest energy particle-hole state, which will in general be smaller than 
the optical gap. This corresponds to putting the particle at the bottom of the single particle band and the hole at 
the top 

Amou = U + 4Jl + ps + ilelik) -/,{k'))^,^^Ck,e',,) ■ (35) 

In Figure 12 we present a caricature of the band structure around n = 4 based on a an explicit non-interacting band 
structure calculated for K4C60 by Gunnarsson et alm^ which is rescaled by a factor 1/3 for the particle band and 
2/3 for the hole band. For this band structure, and any other where the max and min are not at the same k- vector, 
we find an indirect gap where the optical gap is larger then the Mott gap and that the lowest energy particle-hole 
excitations have non-zero momentumiS^ 

r HP r X M 
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FIG. 2: Sketch of the band structure at doping n — 4: (K4C6C^^) showing the indirect gap structure with optical gap Aopt and 
Mott gap Aiviott. The thin dashed lines indicate the centers of the particle and hole bands and the energy scales on the right 
are defined with respect to these. Ao = U + 4Jl + ^Js is the gap in i = limit. 



A word of caution may be appropriate in considering this figure, namely that the single particle states in the bands 
are only well defined within our theory for a small density of such states. We depict the lower "band" as filled with 
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single particle states, but the real entities are only the holes in this band. This is a strongly interacting system and 
the analogy with a weakly interacting semiconductor has limitations. For instance it is quite obviously nonsensical to 
fill up the particle band with a density of more than two particles because that would correspond to a total electron 
density of more than six. We will return to issue of doping away from the Mott insulator in section HVI 

For comparison with other models of the Mott transition on degenerate Hubbard models it is useful to write down 
a more general expression for the Mott gap. If wc make the reasonable assumption that the top and bottom of the 
band structure are roughly the same magnitude W/2, where W is the bandwidth of the tight binding Hamiltonian h, 
we can write 

AMott«t/e//(4)-iM/, (36) 

where /7e//(4) = U + 4Jl + § Js- This expression is in sharp contrast to calculations on iV-band Hubbard models 

without the multiplet sphtting terms, and L^, where forms such as Aiviott ~U — NW or Aiviott ~U — \/NW has 
been suggested.^"* The intuitive motivation for the N dependence is an increase in the kinetic energy of the particle 
hole state due to the additional hopping channels. In this model we see a different behavior since the number of 
hopping channels are limited by the strong spin dependent on-site interactions. 

C. Why is A3 Ceo metallic? 

One of the most striking facts about the fuUerides is that the A3C60 materials are generally metallic given that 
the even integer filling materials are likely large U insulators. The problem of odd- integer filling is significantly more 
complicated than that of even integer filling as presented above. The reason for this is that even in the strong coupling 
limit, Eq. the ground state consists of states (n = 3, L = 1, S* = i) with spin and orbital degeneracies. This problem 
resembles the half-filled single band Hubbard model with a highly degenerate ground state which will be split to order 

/U . The ground state may then have magnetic and/or orbital order. 

However, if we assume that the putative insulating ground state is not ordered so that the hopping of particles, 
(n = 4, L = 0, 5* = 0), and holes, {n — 2, L ~ = 0), is not frustrated by the spin interactions we can derive a Mott 
gap in analogous fashion to that for n = A above which reads 

AMott(n = 3) > C/-4Ji-^J5 + (^£Kfc)-^4(^^')Wfc.,,,) 

« Ueff{i)~\w, (37) 

where W again is the non-interacting bandwidth of h. If the ground state has significant magnetic or orbital corre- 
lations we expect the gap to be bigger because of a lower ground state energy and frustration of the motion of the 
particle and hole. 

Compared to the expression 1361 for the Mott gap at n = 4 we note an increase from 14^/2 to 2W/3 in the kinetic 
energy of the particle and holes due to the larger phase space allowed for hopping. In A3C60 where the bandwidth is 
around .5eV it appears that this difference will not be large enough to close the .5eV Mott gap seen in A4C60. 

However, in addition there is in this model also a more distinct difference between even and odd integer filling, 
namely the sign change of the J5 and terms between the effective Hubbard repulsion Ueff{A) = U + 4Ji + f and 
Ueff{3) = U -4Jl- fJs- This difference comes from the fact that the lowest energy particle and hole excitations 
are (L = 0, S* = 0) at odd integer filling while they are {L = 1, S — ^) at even integer. If the multiplet splitting 
interactions are large enough they could certainly destabilize the Mott insulating ground state at odd integer filling. 
For instance, if Js and Jl are very large such that Ueff{3) < and |J7e//(3)| >> t we would have a spinless Bose 
liquid consisting of an equal number of two and four particle singlets which could only propagate to second order in 
tm^ At intermediate coupling Ueff{i) ^ t we would expect some correlated metallic state with most of the spectral 
weight in low spin configurations of the two, three and four particle states, allowing for hopping to first order in t. 
Evidence for the formation of singlet configurations on short time scales in the metallic fullerides Na2CsC6o, R-BaCgo 
and the quenched cubic CsCgo {n — 1) have been presented from NMR spin-lattice relaxation measurements i^LS^ 

In addition, we know that A3C60 is quite close to a metal insulator transition. It has been found that intercalating 
ammonia into the crystal can cause a transition into an insulating magnetically ordered phase. The main effect 
here is presumably the expansion of the lattice and corresponding decrease in the bandwidth, although the crystal 
symmetry is also reduced which may be important in facilitating a magnetically ordered ground state. An important 
consequence of these findings if interpreted through our model is that for A3C60 C4//(3) > because the ground state 
for small t is magnetic and that |C/e//(3)| ~ t because changes in the magnitude of t can induce a metal insulator 
transition. 
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D. Experiment 

A basic observation from experiments is that there appears to be two distinct energy scales in these materials. 
Probes that are sensitive to spin, in particular NMR, are consistent with a spin gap of around 50-lOOmeV while 
optical conductivity sees a larger charge gap of around SOOmeV. Here we will look at the consistency of the model 
with these observations and make a fit to estimate our microscopic parameters J7, Jl and Js- 



1. Optical Gap 



The charge gap is seen in optical conductivity as a depletion of the low energy weight in K4- and Rb4C6o below 
roughly 500meV."' By inspection of Figure |21 together with the expression 1341 for the optical gap we get the following 
estimate 

Aopt(if4C6o) ~ Ueff{4) - 200meV w SOOmeV', (38) 

where again C/e//(4) = U + 4Jl + ^Js- Solving Eq. EHI gives C/e//(4) ~ 700meV. 

We can compare this to AaCgo where we expect the charge gap to close. Using expression 13 71 for the charge gap at 
n = 3 together with a bare bandwidth of 600meV gives 

AMott (AsCeo) ~ Uef f (3) - 400mey < , (39) 

with Ueff{3) = U — AJl — ^Js- Together with the estimate Ueff{3) > as discussed in Sec. IIII CI we thus find the 
rough estimate of < Ueff{3) < iOOmeV. 

Combining the values for Ueff{^) and Ueff{3) we can now estimate the microscopic parameters of the model. We 
find U < 550meV and 150meV < § Js + AJl < 350meV. 



2. NMR, l/Ti 

Various probes^ have detected a thermally activated magnetic susceptibility in K4C60 and Rb4C6o and more 
recently also in Na2C6(j^. This has been interpreted as evidence for a singlet-triplet gap of Jahn- Teller distorted 
molecules where a a molecule is thermally excited from the JT ground state singlet to the triplet which then acts 
as a local moment. In the model presented here it is natural to assign such experimental signatures of gapped spin 
excitations to the spin modes or magnons which are a necessary part of the spectrum of the non-magnetic Mott 
insulator. 

We will be focusing on measurements of 1/Ti, the spin lattice relaxation rate, deriving the temperature dependence 
of the relaxation by the magnons in the limit T << A^. 

The probability of a transition between nuclear spin states with z-component m' — m ± 1 due to a two magnon 
process which scatters a magnon with i ~ {L^-,S^) and momentum k to i' — {L'^,S'^) and k' is given by Fermi's 
golden rule as 

= Y^{m,n^k,n\,k'\V\m\n^k-l,n\>k' + l)\^S{E,k- E,,k'), (40) 

ik,i' k' 

Here nn, is the magnon number operator, the interaction V = AI ■ Sq with hyperfine coupling A, nuclear spin / and 
electron spin at the nuclear site So. We have dropped the tiny, typically lO^^eV, Zeeman splitting of the nuclear 
spin. 

The electron spin operator will act as a ^/^gz+fm^nj/j^r.L^s^ with some small prefactor given by the overlap of 
these states with the single atom, where bl ^^g^ creates a triplet excitation. Ignoring the details of the matrix element 
between the different triplet states one finds that the general magnon matrix element is given by (nj/fc/ -|- l){nik), 
where {riik) = {e^^^'^ — \)~^ is just the Bose occupation of the magnons. Finally the relaxation rate is given by^ 
_ 2^ mm' — mm'^^ni — .^y]^g]-g ^^iQ maguou part of Wmm' clearly is independent of m and to'. We thus arrive at the 
final expression 

- I de . (41) 



A, 



sinh^(/3e/2) 
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where we have converted the sums to integrals by introducing the density of magnon states N{e). Here is the spin 
gap, i.e. the lower edge of the magnon band, and Wniag is the magnon band width. 

For T << As we can replace the 1/ sinh^(/3e/2) by 'ie~^'^ and the integral is dominated by e ~ A^. Assuming a 
quadratic dispersion at the band edge we get N{e) ^ v^e~~A7 for e > A^. By change of integration variable we arrive 
at the temperature dependence 

l/Ti - r2g-A,/T ^ T << A, (42) 
Figure 131 shows a fit of this model to ^■^C l/Ti data on K4C60 and Na2C60'^ 
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FIG. 3: Fit to ^^C NMR data by Brouet et ali^. The thick lines are fits to l/Ti ~ t^^-^s/t ^j^j^ ^ 24OK for K4C60 and 
As = 710A' for Na2C6o. The thin lines are fits to l/Ti ~ e"'^^'''^ with corresponding = 660A' and 1260if. (The bump 
around 180K for Na2C6o is a presumably a molecular motion peak.) 

We get an excellent fit to the below room temperature activated behavior with a value As — 24QK sa 25meV 
for K4C60- For comparison we also show fits to a model of localized triplet states corresponding to Wmag = and 
1/Ti ~ e"^"/"^. This is the fit used in the the experimental wor k?.»M§. which is based on a model of a static uniform 
Jahn- Teller singlet-triplet gap. In this intermediate temperature regime it is difficult to tell which fit is best, in 
particular considering the fact that both models clearly fail at higher temperatures where the relaxation rapidly 
saturates, and we conclude that the NMR data cannot resolve the two scenarios. 

Nevertheless, using the value A^ — 25meV for the spin gap we may estimate the microscopic parameters. Given a 
bandwidth Wmag which we assume to be symmetric around the center we get 

As = A - iymag/2 , (43) 

where A — 2Jl + 2 J5 is the t — spin gap. (There is an additional corrections to the spin gap of order P/U which 
is a shift of the ground state energy which should be included in a more rigorous treatment.) A very rough estimate 
of the bandwidth may be given by the degeneracy of the spin modes Wmag = 9i^/C/. With U « 500meV as derived 
from charge gap and taking t « WOmeV form band structure calculations gives Wmag ~ 200meV. Collecting into 
Equation 221 for the spin gap gives A = 2Jl + 2Js ~ 125meV which seems in reasonable agreement with the estimate 
IbOmeV < | + 4 J/, < 350meV from the charge gap. This order of magnitude agreement for the coupling constants 
Js and Jl is certainly encouraging in that it comes from experiments on two apparently separate physical quantities. 

A fit to 1/Ti for Na2C6o gives a larger gap of around TOOK. Within our model the t = spin gap A is the same 
for rt = 2 and n = 4 so the differing spin gaps are somewhat unexpected. However, since the crystal structure is 
different, the tight binding Hamiltonians of these materials may be very different and consequently the bandwidth of 
the spin modes. In fact, Na2C6o is fee while K4 and Rb4 are body centered tetragonal. The natural interpretation 
for the variations within this model is thus variations of the magnon bandwidth. Along the same lines we note the 
behavior of 1/Ti in Rb4C6o under pressure where it is found that the activated behavior is replaced or coexists with a 
non-activated component related to gapless excitations.^ It has been suggested that this is related to a closing of the 
Mott gap due to the expected pressure induced increase of the bare bandwidth. Within our model we find a possible 
alternative interpretation in terms of a closing of the spin gap. 

Above room temperature the activated behavior stops and the the relaxation rate saturates. Within our simple 
non-interacting model for the spin modes we cannot expect to be able to address the high-temperature behavior when 
a significant number of modes are excited. A more sophisticated treatment requires us to properly account for the 
interactions between the spin modes as well as the exclusion statistics that are ignored in the single particle picture. 
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The saturation could also be related to molecular degrees of freedom at higher temperatures which are completely 
neglected in our modeL'^" However, we note that this rapid saturation is very reminiscent of the behavior of l/Ti in 
spin ladder materials with gapped magnons^, where it is believed to have a purely electronic origin'^''. 

IV. THE DOPED MOTT INSULATOR 

We will now look at the problem of an incommensurate particle density away from the Mott insulators at even 
integer filling n = 2 and n = 4. This is obviously a much more difficult task because in the strong coupling limit the 
ground state will be highly degenerate. As a concrete example, at a doping n = 2 + x {x < 1) the ground state at t = 
is the set of states with x 3-particle states (3, 1, ^) and (1 — x) 2-particle states (2, 0, 0) at arbitrary positions in space. 
Introducing t by means of the the effective Hamiltonian (|17(l we find to first order in t nearest neighbor interchange 
of the 2- and 3-particle states and to second order in t spin and orbital exchange terms between nearest neighbor 
3-particle states. This can be described by a generalized t — J-model including a no double occupancy constraint 
because only 2 and 3 particle states are allowed. 

At first glance this may appear to be an even more difficult problem than that of a doped antiferromagnet because 
of the additional orbital degrees of freedom. However, in the low density limit, x << 1, it is in fact considerably 
simpler than the doped antiferromagnet because only the doped particles (or holes) have internal spin and orbital 
degrees of freedom. In the doped antiferromagnet the scenario is just the opposite with a large number 1 — x of 
spinful particles and a small number x of spinless holes. This of course gives rise to the very complex behavior in such 
systems where the spin interactions J can compete with the hopping t even in the limit J ^ /U « t because the 
important energetics is given roughly by xt and (1 — x) J. Here for the doped non-magnetic Mott insulator a similar 
consideration would lead us to compare xJ with xt because it is the doped particles or holes that carry both the spin 
and momentum. Effectively we are thus looking at the low density (heavily doped) limit of a t-J model. 

We will completely neglect the nearest neighbor exchange interactions as well as the no double occupancy constraint 
and only consider the single particle physics. Certainly, for the problem of a single particle or hole doped into the 
non-magnetic Mott insulator this is completely rigorous and again in sharp contrast to the problem of a single hole 
in an antiferromagnet where interactions obviously cannot be neglected. Even this single particle physics has some 
interesting implications for the doping dependence in the metallic fuUerides. 

A. Small Fermi surface 

The problem of a single particle or hole was addressed already in Section IIII B II in connection with particle- hole 
excitation. There we showed that the single particle or hole spectra are equivalent to the non-interacting Hi = 
spectrum up to a rescaling by a factor 1/3 or 2/3. At least in the very low density limit, n~2 + xoTn = 4: + x with 
|a;| << 1, we expect these single particle states to give a qualitatively accurate picture by filling x such states. 

In particular this implies a "small Fermi surface" where the number of delocalized charge carriers is proportional 
to the number of doped holes or particles x and not the total filling n. The remaining degrees of freedom are frozen 
below the Mott gap. One important consequence is that the density of states at the Fermi surface for small doping x 
will be given by the density of states at the band edges of the of non-interacting problem by the simple relation 
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A detailed picture of what happens at larger doping x — > 1 as we approach odd integer filling is beyond our methods. 
However, a naive extrapolation of the results valid for small x all the way to x = 1 gives a density of states as shown 
in Figure 01 where the density of states is generally peaked at odd integer filling as a consequence of the rapid decay 
toward the effective band edges at even integer filling. For reasons discussed in Section IIII CI we have to be in an 
intermediate coupling regime where the system is metallic at odd integer filling for this extrapolation to have any 
credibility. 

There is an interesting experiment that corroborates the small Fermi surface picture in Ceo which is the variation 
of density of states in Na2CS:cC6o (0 < a; < 1) corresponding to a doping range 2 < n < 3i It was estimated from the 
Pauli susceptibility that for samples with n = 2.25, 2.5, 2.75, 3 the density of states varies as 5, 7, 11, 15 eV~^ (both 
spins). Correspondingly, Tc drops rapidly from 12K at n = 3 to 7K at n = 2.75 and to < Q.'bK at n = 2.5. Certainly, 
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this behavior seems consistent with the scenario sketched in figure^ where the density of states drops rapidly as the 
effective Hubbard band edges are approached at even integer fiUing. 

The failure of naive band theory in the presence of strong local repulsion is a consequence of the large inherent 
charge fluctuations of such an uncorrelated state of delocalized electrons. It may be illuminating to recall some 
of Hubbard's original work on the topic of narrow band systemsi^ For an m-fold degenerate band at filling n the 
probability PAr(m, n) of having N electrons on a particular atom (molecule) is given by 



P^Kn)=(™ (-^(1--)-^, (45) 
\i\ J m m 

where (^) is the multiplicity of atomic states with N particles. The rms fluctuation is given by (AAr)RMS = 
•y/n(l — n/m) which has a maximum (A7v)rms — \/rn/2 at half filling n — m/2. Clearly, at finite doping there 
are significant charge fluctuations of order one which cost an energy of order U per site in an uncorrelated state and 
which grows with the degeneracy. We can get an estimate of the energy cost of the charge fluctuations for our model 
by comparing the potential energy E = {Hj) in an uncorrelated state which is the ground state of the kinetic energy 
h with that given by the small Fermi surface state which is the ground state of the potential energy Hi. The latter 
is at filling n given by pjv N-particle states and pn+i N + 1-particle states in the lowest energy multiplet where 
< n < iV + 1 and pnN + pn+i{N + 1) = ti. The potential energy of this state is given by 

{Hi),,,,=PnEo[N)+pn+iEo{N +1) , (46) 

where Eo(N) is the energy of the ground state multiplet with N particles. The expression 1)45(1 turns into 

PN,LAn) - (2L + 1){2S + 1)(^)^(1 - ^)«-^ (47) 

for this model where the A^-particle multiplets are split according to L and S and the corresponding potential energy 
for the uncorrelated state is 

(i^/)u„corr= PN,L,sin)E{N,L,S). (48) 

N,L,S 

FigureElshows AE = (i?/)uncorr — (^^/)corr as a function of doping and in units of Ueff{2) = t/gjy (4) = U+4:Jl + 3/2Js 
for two cases Jl = Js = and Jl = -lU, Js = -211. 
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FIG. 5: Potential energy difference, in units of Ueff{2) = U + AJl + 3/2 Jg, between an uncorrelated and correlated ground 
state as a function of filling n. The solid line is for Js = 2Jl ~ -lU and the dashed line forj^ — Js = 0- 

We see that in a wide doping range around half fiUing there is a significant energy cost due to charge fluctuations 
in an uncorrelated state. Given that this leads to the Mott insulating behavior at even integer filling we would also 
expect significant correlation effects in the metallic regions 2 < n < 4, consistent with the small Fermi surface scenario. 

Hubbard's treatment of narrow band systems used a real space Greens function method which is exact in the zero 
bandwidth t = case but which depend critically on neglecting correlations between electrons on different atoms in 
the finite bandwidth case. He found that in the narrow bandwidth limit the original non-interacting band splits into 
a large number of bands which correspond to transitions between states with particle number differing by one. In 
addition these Hubbard bands reflect the non-interacting band with density of states which are some functional of the 
non-interacting density of states. Our result for the doped Mott insulator is essentially a special case of Hubbard's 
results for which we can solve for the quasi particle spectrum exactly in the low density limit where we can neglect 
interactions between the doped particles or holes . 

B. High Temperatures 

In the previous section we discussed the ground state properties of the doped Mott insulator close to even integer 
filling. We found a band of single particle states which up to a rescaled hopping are equivalent to the states of 
the non-interacting problem. At low temperatures we thus expect a simple metallic behavior with band like charge 
transport. Here we will speculate on some of the interesting physics which could emerge from the model at higher 
temperatures. 

The derivation of the single particle Hamiltonians, Eq. 1261 and [771 for the particle and hole states depend crucially 
on the fact the Mott insulating ground state is of the simple non-degenerate form \GS, 4/2) = Y[r I"- ~ 4/2, L = 0, S ~ 
>r and that we can ignore the higher orbital and spin multiplets in the 2 and 4 particle molecular spectra. At elevated 
temperatures of order the spin gap this assumption is no longer justified as the the spin and angular momentum 
modes discussed in Section llll Al are thermally occupied and we should consider their effect on the quasiparticle states. 
(The activated behavior discussed in Sec. IIIID 21 which has been linked to gapped triplet excitations is seen in NMR 
(1/Ti) also in the metallic Na2CsC6o and RbaCgoi^^) These, in fact, interact very strongly with the quasiparticles 
in a quite non-trivial fashion. In the presence of such an excited spin or orbital state the nearest neighbor hopping 
integrals described by Equation 1221 turn into some more complicated expression given by for instance 

(4, 1, 1, L'\S''\r{3, I', s'\r>{J2 t73clj,aCr',r,a)\3, I, s)r\i, 1, 1, i^ S')r' , (49) 

jj', Cl- 
in the case of a hole hopping to a site occupied by a spin triplet. Depending on the configurations of the spin triplet 
states the hopping of the hole (|3, ^,s)) may be completely suppressed or it may require a spin flip. Certainly the 
nearest- neighbor hopping integrals would be completely altered from the simple form of Equation 1221 

It is not obvious how to model this problem but the most naive scenario might be to ignore the hopping of the 
spin and orbital states and replace them by thermally excited impurities causing disorder in the hopping of the 
particles or holes. Such a model would be similar to that suggested by Varma to explain the paramagnetic insulator 
to ferromagnetic metal transition in lanthanum manganites (giant magnetoresistive compounds). '^^ 

One may speculate that this temperature activated off-diagonal disorder could destroy the band like motion of 
the charge carriers and possibly be related to anomalous properties of A^Cqq at elevated temperatures such as the 
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evidence for localization from NMRSi, the disappearance of the Fermi edga^ and the non-saturation of resistivity 
and the correspondingly extremely short mean free paths^. 

V. CONCLUSIONS 

We have studied an orbitally degenerate 3-band Hubbard model with additional multiplet splitting on-site inter- 
actions JsS^ and JlL^ which favor low spin and low orbital angular momentum. We use the effective Hamiltonian 
method in the strong coupling limit U » Js » Jl >> t perturbatively to second order in t. At even integer filling, 
n = 2 or n = 4, this model is a insulator with a non-degenerate ground state where the electrons at each site occupy 
the L = and S = configurations and with distinct spin and charge gaps. The trivial ground state allows for a 
simple single particle description of spin and charge excitations. The lowest energy spinful excitations is a band of 
magnons with a bandwidth Wmag ^ /U and a gap As = 2 Jl -|- 2 Js — 0{t^ /U). A single particle or hole doped into 
the Mott insulator are described by the non-interacting tight-binding Hamiltonian but with an overall rescaling by a 
factor 1/3 or 2/3 of the hopping integrals and corresponding bandwidth. The latter allows for a detailed description 
of the particle-hole excitations and the corresponding charge gap is given by Ajviott ~ U -\- 4Jl -|- 3/2 J5 — W/2 in 
terms of the band width W of the non-interacting Hamiltonian. 

Close to the Mott insulator, at filling 2 + x or 4 -|- a; with |a;| << 1, we find a metallic state with a "small Fermi 
surface" where the density of charge carriers is given by |a:| and a density of states which is simply a renormalization 
by a factor 3 or 3/2 of the density of states at the band edges of the non-interacting band structure. Consequently, 
in three dimensions the density of states will in general increase rapidly with the the doping x. 

In this model there is also a distinct difference between even and odd integer filling which follows from the simple 
fact that an odd number of electrons cannot form a spin singlet. From this follows that the effective on-site repulsion 
is given by [/eff(n = 2/4) + AJl + 3/2Js at even integer filling and by C/off(n = 1/3/5) = U - AJl - 3/2 Js at 
odd. Consequently, depending on the magnitude of Jl and Js, the Mott gap may be significantly reduced or vanish 
at odd integer filling. 

The properties of this model are strikingly similar to the phenomenology of the fuUerides A„C6o with 2 < n < 4. The 
non-magnetic Mott insulator at even integer filling with a small spin gap and a larger charge gap, the even/odd effect 
at integer doping where A3C60 is generally metallic, as well as the rapid suppression of the DOS and the corresponding 
superconducting transition temperatures as the filling approaches even integer. We do a fit of the model to the charge 
gap from optical conductivity and the spin gap from NMR l/Ti in K4C60 which appear consistent with values of Js 
and Jl of around 50 — lOOmeV. 

There is a number of interesting open questions about the model and the possible implications to alkali doped Cgo- 
In particular we need a better understanding of the physics at odd integer filling on the metallic side of the Mott 
transition. Can this state have a superconducting ground state even though, as evidence suggest, [/e//(3) > such 
that there is no bare attraction in the way envisioned by Chakravarty and coworkers as an electronic mechanism 
of superconductivity? In fact, also at even integer filling there have been intriguing suggestions of an intermediate 
superconducting state in the metal-insulator transitionpiSiii Another issue is the properties of the model at elevated 
temperatures approaching the spin gap where we have found that the spin and orbital modes interact strongly with 
the charge carries and may significantly effect the simple band like charge transport. 

APPENDIX: SECOND ORDER CANONICAL TRANSFORMATION. 

Here we derive the expressions for the second order terms in the effective Hamiltonian Eq. 1171 Starting with Eq. 

cni 

where Si is given by 

- 7 l,rr 

(A.1) 

Our purpose is to construct ^2 such that it cancels terms which connect different energy subsectors from the 
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commutators /i^^"'^]. We have 

il.rr' 

{ry).,{r".r"') ^ ^ ^ ^ 

J l,rr 

- E j^(^)'!!j^(„) ^c^'^'''''"[l«)>-^'(H.rs|c).v--(rf|.v»], (A.2) 

{T.T\r") 

where we have used the fact the operators commute if there is no overlap between sites. Using the complete set of 
three site states 1 = \a) „> r" (p\rr' r" for the three site interaction we arrive at 



CL 1 rr' L l:'"^' \ 

We now split this into a part SE{a)E{b) which is diagonal in energy and (1 — (5£;(Q)£;(f,)) which connects different 
energy sectors. The latter part we can cancel by solving for ^2, schematically 

[S2,Hj] = + /ll])off-diagonal ■ (A.4) 

Using an ansatz ^2 = J2{r,r')^2'^lbW)rr'{b\rr' + J2{r,r'.r")^2'^lb"W)rr'r"{b\rr'r" and using the fact that 
[\a)rr'r" {b\rr'r" , Hj] = {E{b) — E{a))\a)rr'r" {b\rr'r" , wc Can solve for ^2 ~ 0{P / AE). Apart from the cancellation, 
S2 will only contribute to higher orders. 

Having done the cancellation we are left with terms that are diagonal in energy. Clearly [Si , does not contribute 
to this because h'^ is diagonal in energy and S ^ is strictly off-diagonal. The result for the remaining second order 
terms is 

/l^] diagonal = ^9 E ~ TP(c) - EU) ~ ^ E{a]E{b)\a) rr' {b\rr' 

(r,r'),c ^ ' ^ ' 

- 2 E ~ £;(c) - E{a) ^^^l^i'>')E{b)W)rr'r"{b\rr'r" , (A.5) 

{r,r' ,r") ,c 

which simplifies to the final expression given in ll7l 
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